Collective states of interacting Fibonacci anyons 
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We show that chains of interacting Fibonacci anyons can support a wide variety of collective ground states 
ranging from extended critical, gapless phases to gapped phases with ground-state degeneracy and quasiparticle 
excitations. In particular, we generalize the Majumdar-Ghosh Hamiltonian to anyonic degrees of freedom by 
extending recently studied pairwise anyonic interactions to three-anyon exchanges. The energetic competition 
between two- and three-anyon interactions leads to a rich phase diagram that harbors multiple critical and gapped 
phases. For the critical phases and their higher symmetry endpoints we numerically establish descriptions in 
terms of two-dimensional conformal field theories. A topological symmetry protects the critical phases and 
determines the nature of gapped phases. 

PACS numbers: 05.30.Pr, 73.43.Lp, 03.65.Vf 



Two-dimensional topological quantum liquids such as the 
fractional quantum Hall (FQH) states harbor exotic quasipar- 
ticle excitations which due to their unusual exchange statis- 
tics are referred to as anyons. Interchanging two anyons may 
result in not only a fractional exchange phase, but may also 
give rise to a unitary rotation of the original wave function in 
a degenerate ground-state manifold. This latter case of non- 
Abelian statistics is proposed to be exploited in the context 
of topological quantum computation |Q1 0]. Intense experi- 
mental efforts |3[] are currently under way to demonstrate the 
non-Abelian character of quasiparticle excitations in certain 
FQH states, as proposed theoretically HI lit]. 

Given a set of several non-Abelian anyons we can ask what 
kind of collective states are formed if these anyons are inter- 
acting with each other. A first step in this direction has re- 
cently been taken by studying chains of "Fibonacci anyons" 
with nearest-neighbor interactions [6]. Fibonacci anyons rep- 
resent the non-Abelian part of the quasiparticle statistics in 
the k = 3 Zfc-parafermion 'Read-Rezayi' state J3l, an effec- 
tive theory for FQH liquids at filling fraction v = 12/5 [7]. 
A single Fibonacci anyon carries a topological charge r. Two 
such anyons may combine ("fuse") so the pair has charge r or 
has no charge, which is denoted 1. This is analogous to two 
St/ (2) spin- 1/2 's combining to either spin-1 or spin-zero total 
spin. A two-anyon interaction assigns different energy to the 
two possible charges of the pair, just as a Heisenberg exchange 
interaction does for the two possible total values of spin of a 
pair of SU(2) spin-l/2's. For a chain of Fibonacci anyons 
with a uniform pairwise nearest-neighbor interaction of either 
sign it has been explicitly shown [6] that the Hamiltonian has a 
topological symmetry, which was predicted to stabilize one of 
the gapless phases. In this manuscript, we give a broader per- 
spective on possible collective phases of interacting Fibonacci 
anyons and phase transitions between them. 

The specific model we will focus on has, in addition to the 



two-anyon term, an additional three-anyon interaction both of 
which may arise from tunneling 0]. We find a rich ground 
state phase diagram that harbors multiple critical, gapless 
and gapped phases. We also mention in passing other one- 
dimensional models that we have investigated and will report 

J— 1 

on in a forthcoming longer paper |9|]. The topological sym- 
metry, introduced in Ref. 6, that measures the topological 
flux through a ring of Fibonacci anyons plays an essential role 
in determining the nature of the observed phases and phase 
transitions. In particular, we find that this topological sym- 
metry protects all the critical phases against spatially uniform 
local perturbations. These extended critical phases can be de- 
scribed in terms of 2D conformal field theories (CFT) with 
central charges c = 7/10 and c — 4/5 and can be respectively 
mapped exactly onto the tricritical Ising and 3-state Potts crit- 
ical points of the generalized hard hexagon model lallOIl- At 
the phase transitions out of the tricritical Ising phase into ad- 
jacent gapped phases the system exhibits even higher symme- 
tries which we identify as fefracritical Ising and 3-state Potts 
critical points. Remarkably, this demonstrates that these 2D 
classical models share an identical non-local symmetry which 
is the classical analog of the topological symmetry in the ID 
quantum chains. At the transition into the gapped phases this 
topological symmetry is spontaneously broken which results 
in a non-trivial ground-state degeneracy in the gapped phases. 





FIG. 1 : (color online) Illustration of a) the alternating chain and b) 
the generalized Majumdar-Ghosh chain with a three-anyon interac- 
tion term. Shaded enclosures indicate the fusion products that are 
energetically biased by the Hamiltonian. 
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Two-anyon interactions- For a uniform chain of Fi- 
bonacci anyons the Hamiltonian introduced in Ref. |g ener- 
getically favors one or the other of the possible fusion prod- 
ucts of two neighboring r-particles which, by the fusion rule 
t x t = 1 + t, can be either a 1 or a r. The energy of the for- 
mer is lower for a coupling that is termed 'antiferromagnetic' 
(AFM), in analogy to the familiar SU(2) spins-chains, while 
that of the latter is lower with a coupling termed 'ferromag- 
netic' (FM). The underlying Hilbert space is spanned by an 
orthonormal basis of states, each state corresponding to one 
possible labeling of the chain [6] of repeated fusions with r. 
Each site along this chain of fusions has either a 1 or a r, with 
a constraint forbidding two adjacent l's. 

By performing a sequence of local basis transformations 
and projection onto one of the two fusion channels for each 
pair of neighboring anyons, the resulting two-anyon interac- 
tion Hamiltonian can be written as a sum of local 3-site oper- 
ators H = J-2 Y^i H\ which take the explicit form (%' denotes 
the first in a triple of adjacent sites) 

H\ = -Vlrl ~ <T 2 7>TlT - f^TTT (1) 

-</}- 3/2 (\t1t) (ttt\ + hue.) , 

where V a projects onto the state \a), e.g. Viri = |1t1) (1t1| 
and <j> = (1 + y/E) /2 is the golden ratio @|. 

Here we want to explore a larger space of models than that 
given by this uniform chain with only nearest-neighbor two- 
anyon interactions. One way is to let the strength J 2 of the 
interaction alternate along the chain, as illustrated in Fig.QJ). 
Two chains can be coupled to form a two-leg ladder. Another 
way is to add a spatially-uniform three-anyon interaction, as 
indicated in Fig.QJ)), which because of its rich phase diagram 
we discuss in the following. 

Majumdar- Ghosh (MG) chain- Three SU(2) spin-l/2's 
can combine to a total spin 3/2 or 1/2. For a uniform SU (2) 
spin- 1/2 chain, Majumdar and Ghosh showed that an AFM 
coupling (favoring total spin 1/2) for each set of three neigh- 
boring spins gives rise to a gapped phase with the two possi- 
ble dimer coverings being the exact ground states [11]. In the 
same spirit we have asked what possible phases can be stabi- 
lized by a spatially uniform three-particle interaction term in 
our 'anyonic generalization of the SU(2) Heisenberg model', 
i.e. a term that energetically favors each set of three adjacent 
Fibonacci anyons to fuse together ]8|] into either a 1 or a t, as 
illustrated in Fig.QJ)). Like the pairwise interaction term such 
a three-particle interaction term respects both the translational 
and topological symmetries. We find that the energetic com- 
petition between such two- and three-anyon interactions gives 
rise to the rich ground-state phase diagram shown in Fig. [2] 
which we discuss in some detail in the following. Similar to 
the derivation for the pairwise interaction term <j2j we can ob- 
tain a local form H = J3 £\ H3 °f me three-anyon interac- 
tion term by a sequence of basis transformations and projec- 
tions which then takes the explicit form of a 4-site interaction 
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FIG. 2: (color online) The phase diagram of our anyonic Majumdar- 
Ghosh chain on the circle parameterized by 9, with pairwise fusion 
term J2 = cos# and three-particle fusion term J3 = sin#. Besides 
extended critical phases around the exactly solvable points (9 = 0, n) 
that can be mapped to the tricritical Ising model and the 3-state Potts 
model, there are two gapped phases (grey filled). The phase tran- 
sitions (red circles) out of the tricritical Ising phase exhibit higher 
symmetries and are both described by CFTs with central charge 4/5. 
In the gapped phases exact ground states are known at the positions 
marked by the stars. In the lower left quadrant a small sliver of an 
incommensurate phase occurs and a phase which has ^-symmetry. 
These latter two phases also appear to be critical. 

between consecutive labels along the chain of fusions 

H\ = V T \t\ + V\ T l T + V T TTl + VlTTT + 2(j) 2 V TT TT + 

0" 1 (V TlTT + V T rlr) - </>~ 2 (\tt1t) (t!tt\ + h.C.) + 

0~ 5/2 (\t1tt) (tttt\ + \tt1t) (tttt\ +h.c.) , (2) 

where the site i denotes the first position in each 'quad' of 
sites. The full Hamiltonian with competing fusion terms then 
becomes Hj 2t j 3 = 2^1(^2-^2 + J^H^), where we parameter- 
ize the couplings by the angle 9 as J2 — cos 9 and J3 = sin 6. 
We study periodic chains of L anyons. 

The phase diagram of this model, shown in Fig.|2j exhibits 
two critical phases that contain the two exactly solvable points 
(9 = 0, 7r). These extended critical phases can be described 
by 2D conformal field theories and are thereby related to 2D 
classical critical points to which an exact mapping was estab- 
lished at the two solvable points [6;]. For AFM pair interaction 
( J2 > 0) this is the tricritical Ising model (c = 7/10), while 
for FM pair interaction ( J2 < 0) it is the critical point of the 
3-state Potts model (c = 4/5). In particular, we note that 
the critical phases found at the exactly solvable points are sta- 
ble upon introducing a small three-anyon fusion term. While 
the J3-term respects both translational and topological sym- 
metries, all translational invariant operators with scaling di- 
mension < 2 at the exactly solvable points are found to break 
the topological symmetry [6]. This shows that the topologi- 
cal symmetry protects the gaplessness in the vicinity of these 
points, somewhat analogously to the much-discussed notion 
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FIG. 3: (color online) Energy spectrum at the S^-symmetric point 
{9 = 0.1767r) of the Majumdar-Ghosh chain. The energies have 
been shifted and rescaled so that the two lowest eigenvalues match 
the CFT scaling dimensions. The open boxes indicate the positions 
of the primary fields of the parafermionic subset of the c = 4/5 CFT; 
fields with topological flux r are marked. The open circles give the 
positions of multiple descendant fields as indicated. 

that a topological symmetry protects a ground-state degener- 
acy in a gapped topological phase in 2D space. 

For large three-particle interactions these critical phases 
eventually give way to other phases, such as the two distinct 
gapped phases indicated by the grey shaded arcs in the phase 
diagram. Remarkably, the transition to the gapped phase 
from the tricritical Ising phase when both interaction terms 
are AFM ( J2 > and J3 > 0) apparently has an "emergent" 
S3 (3-state Potts) symmetry. Our numerical analysis shows 
that this transition occurs at 8 = 0.1767T and is described by 
the parafermion CFT with central charge c = 4/5, indicative 
of an additional S'3-symmetry at this point. Fig. [3] shows the 
rescaled energy spectrum at this critical point whose (univer- 
sal) low-energy part is in spectacular agreement with the CFT 
predictions. Note the relevant operator with zero momentum, 
zero flux and scaling dimension 4/3, which breaks the S3- 
symmetry. It is the leading operator present in the Hamilto- 
nian away from this special point, and drives the system into 
either the gapped or the tricritical Ising phase. In the gapped 
phase, the topological symmetry is spontaneously broken and 
the resulting ground state, which has zero total momentum, is 
two-fold degenerate in the thermodynamic limit. In the tricrit- 
ical Ising phase, the Z2 sublattice-symmetry breaking order 
parameter corresponds to a more relevant continuum operator 
than the topological order parameter, while it is the state cor- 
responding to the Z2 order parameter which acquires a higher 
energy in the gapped phase, where only the topological sym- 
metry is broken. At the transition, both order parameters are 
degenerate, see Fig. [3] and together they form the order pa- 
rameter of a critical 3-state Potts model with S'3-symmetry. 

In the case of FM three-particle interaction J3 < the tran- 
sition at 9 = — 0.4727T between tricritical Ising and gapped 
phases is described in terms of CFT by the full c = 4/5 
minimal model representing IU2I1 the fefracritical Ising model. 



FIG. 4: (color online) Energy spectra at the generalized Majumdar- 
Ghosh point with ferromagnetic three-anyon interaction term only 
(8 = 3n/2). With both translational and topological symmetry spon- 
taneously broken there are four ground states for even length chains 
(left panel) with two states each at momenta K = and K = n. 
The gapped low-energy spectrum is a continuum of two-spinon scat- 
tering states (grey shaded). For odd system size (right panel) a dis- 
tinct single-spinon dispersion emerges below a gapped continuum of 
3-spinon scattering states. 

Again we have unambiguously identified the CFT description 
of this critical endpoint by assigning the low-energy states in 
the energy spectrum similar to Fig. [3] (not depicted here). In 
particular, the topological symmetry forces the system onto 
the integrable renormalization group trajectory 111 311 into the 
gapped phase or tricritical Ising fixed point, driven again by 
the relevant operator with dimension 4/3. At this tetracriti- 
cal Ising transition into the gapped phase the system sponta- 
neously breaks both the translational and topological symme- 
tries. As a consequence, we observe a four-fold ground-state 
degeneracy throughout this gapped phase for chains with even 
length. The nature of this gapped phase is best characterized 
at the point 6 = 3ir/2 ( J3 = —1, J 2 — 0) that constitutes the 
anyonic analog of the Majumdar-Ghoshpoint of the spin- 1/2 
Heisenberg chain. At this point the four ground states for even 
L take the exact form 

|-0no-flux) = \T X TT X TT X T . . .) + _1 |t1t1t1 . . .) ± 

\tt x tt x tt x ...)+ 0- 1 |lrlrlr.. .) (3) 
IVv-flux) = \t x tt x tt x t . . .) - \t1t1t1 . . .) ± 

(j)- 1 \tt x tt x tt x ...)- \1t1t1t ...) , (4) 

where t x = <\r x |1) + |t) denotes a normalized su- 

perposition of the states |1) and |r) on a single site. Note 
these ground states have total momenta K = or K = ir, 
indicating the two-sublattice ordering. There are two states 
at each momentum, one with a r-flux and the other without. 
Of course, we can instead make the simpler linear combi- 
nations of these ground states that explicitly break both the 
topological and sublattice symmetries: these four states are 
\t x tt x tt x t . , .), |rlrlrl . . .) and the equivalent states under 
translation by one site. Note that the density of l's [which for 
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FIG. 5: (color online) Schematic energy spectra in the gapped phase 
for 0.19 < 9/tt < 0.30. Below a continuum of scattering states 
(grey shaded) a distinct quasiparticle band forms. The plot combines 
data with 24 < L < 36. The solid lines are a guide to the eye. 



these states is l/(2<p 2 ) and 1 /2, respectively] is a simple order 
parameter that reflects the topological symmetry breaking. A 
pair of these latter states also becomes the exact ground states 
if the Hamiltonian explicitly breaks the sublattice symmetry 
by alternating the two-anyon interaction along the chain as 
shown in Fig. QJ). For a uniform chain in the tri critical Ising 
phase this alternation corresponds to a relevant operator that 
immediately drives the system into a gapped phase. On the 
other hand, the uniform chain in the Potts critical phase has 
Zs sublattice symmetry and no relevant operators with spatial 
period two, so this chain remains critical when a small alter- 
nation is added |@1 • 

The low-energy excitations in the gapped phase around 
the MG point are domain walls between the two sublattice- 
ordered ground states with a low density of l's - similar to 
spinon states in the spin- 1/2 MG chain. For chains of odd 
L the periodic boundary conditions force the presence of an 
odd number of such spinons/walls, and indeed we observe a 
distinct single-spinon dispersion relation over part of the Bril- 
louin zone, as shown in the right panel of Fig. [4] Explicitly 
evaluating the topological symmetry operator [6] we find that 
these one-spinon states all have a r-flux. For even length rings 
we observe a continuum of two-spinon scattering states in the 
spectrum shown in the left panel of Fig. [4] 

At the 5*3 -symmetric transition out of the tricritical Ising 
phase (8 = 0.1767r) the system spontaneously breaks only 
the topological symmetry resulting in a two-fold degenerate 
ground state at total momentum K = in the gapped phase. 



With increasing coupling J3, away from the transition, an en- 
ergy gap opens in the spectrum, with a distinct quasiparticle 
dispersion forming below a continuum of scattering states as 
shown in Fig. [5] Explicitly evaluating the topological sym- 
metry operator [6] we find that these quasiparticle states all 
have a r-flux. For small J3, the lowest energy quasiparticle 
remains at momentum K = n, with both the gap and the 
mass increasing with increasing J3. Near = 0.24-7T the gap 
reaches a maximum value and the mass diverges. The min- 
imum of the quasiparticle dispersion bifurcates and continu- 
ously moves away from K = it, see the lower panel in Fig. [5] 
Eventually, the two minima approach the commensurate mo- 
menta K = 2-7r/3, 47r/3 and as these modes soften the system 
enters the 3-state Potts critical phase at 8 = 0.3167T. Sim- 
ilar to the MG point we can identify a point in this gapped 
phase at which we can determine the exact form of the two 
ground states. For J3 = -| J2, we note that the off-diagonal 
term in Hamiltonian is exactly cancelled by its counterpart 
in the J3-Hamiltonian (|2]), and at this angle, tan (9 = <fi/2, the 
ground states take the explicit form 

\T x f x T x T x ...) + (-1) L _1 \tttt 



IV>n 



-flux/ 



-flux/ 



^ I T~x 7~x 7~x T x ■ 



■)-(-ifh 



(5) 
■(6) 



where f x = (j) 1 / 2 |1) + (f)^ 1 |r) is a single site superposition. 
(A projector onto the states in the Hilbert space is implicitly 
assumed.) Again note that the density of l's is a simple order 
parameter for topological symmetry breaking. 

Finally, we note that when both interaction terms are FM 
the critical 3-state Potts phase gives way to a small sliver of 
an incommensurate phase and then a phase with Z± sublattice 
symmetry. All of these phases appear to be critical or nearly 
critical. In the incommensurate phase, correlations in the local 
density of l's oscillate with a spatial period varying between 
3 and 4 lattice spacings. 

Breaking the topological symmetry.- To further elucidate 
the role of the topological symmetry in determining the nature 
of the observed phases, we can explicitly break this symmetry 
in a modified Hamiltonian. Since the number of l's, Ni, is 
a simple order parameter for this symmetry, we add a term 
hN\ to the Hamiltonian, where h is the field that breaks the 
symmetry. Within the gapped phases the ground-state degen- 
eracy is immediately lifted by h 7^ 0, with the state(s) with 
lower (higher) density of l's favored by positive (negative) h. 
At the point tan 8 = </>/2 these states are precisely the "vac- 
uum" state of l's \tttt . . .) (N% — 0) for h > and the 
state |f2;f x f x f x . . .) (Ni > 0) for h < 0. The whole gapped 
phase for AFM J2 and J3 can thus be identified as a first-order 
transition with a 'liquid-gas' coexistence of l's at h = 0. A 
similar picture, but with broken sublattice symmetry, holds for 
the gapped phase around the MG point. Varying h in the crit- 
ical phases of our phase diagram, there is a weaker feature in 
the ground-state energy E(h): the 'topological susceptibility' 
Xt{h) oc d 2 E/dh 2 diverges with system size at h = 0. This 
indication of a continuous phase transition on varying h con- 
firms that the topological symmetry is a full participant in the 
critical behavior at h = 0. 
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In conclusion, the anyonic quantum chains possess a topo- 
logical symmetry that forces their corresponding 2D classical 
models onto a highly fine-tuned submanifold of their respec- 
tive phase diagrams. The competition of anyonic exchange 
interactions allows one to move within this manifold contain- 
ing a plethora of both, (multi)critical phases such as tricriti- 
cal Ising, 3-state Potts, tetracritical Ising, and various gapped 
phases with spontaneously broken topological symmetry in- 
cluding an analog of the Majumdar-Ghosh point. 
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